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Abstract 

The scaling behavior of fluctuations in cluster size is studied in g = 5 and 7 state Potts 
models. This quantity exhibits scaling behavior on small lattices where the scaling of local 
operators like energy fluctuations and Binder cumulant can not be expected. 
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Rigorous theory of finite size scaling, for first-order phase transitions, put forward 
by Borg et. ah [|l], gave explicit form for the size dependence of various operators. 
According to this theory, for periodic boundary conditions, the operators such as 
specihc heat and Binder cumulant at the transition temperature can be represented 
by a polynomial in 1/L'^ where L is the linear lattice size and d is the dimension. If 
the lattice sizes are much larger than the correlation length, the contribution of the 
higher order terms are negligible 0. The difficulty arises when the correlation length 
is of the order of the lattice size or larger. In this case, higher order corrections are 
necessary and deciding the order of the transition becomes difficult. Even when the 
large lattices are used, higher order terms may create difficulties during the fitting 
process to the simulation data. Such difficulties may be reduced by choosing the 
operators for which the correction terms play less important role. A good example 
for such an operator is the average energy measured at the infinite lattice transition 
point. This quantity has exponentially small correction term which enables one to 
determine the inhnite lattice critical point with great accuracy 

In an earlier work in on the two-dimensional Potts model, it is observed that 
compared to the local operators such as energy and order parameter, cluster related 
(global) operators are more sensitive to the changes in the system going through a 
phase transition. Particularly, the average cluster size distribution may give better 
indication on the order of the phase transition at smaller sizes than that for the en- 
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ergy distribution. The most trustworthy method to test this observation is studying 
the hnite size scaling behavior of the cluster related operators. In this work, the 
hnite size behavior of the cluster size fluctuations in g-state Potts model under the 
light of the hnite size scaling theory, developed for hrst-order phase transitions in 
the past decade [1,2,7,8] will be examined. 


The two-dimensional g-state Potts model is known to undergo a hrst-order phase 
transition for g > 4 H. While the correlation length at this transition, is about 


10 lattice sites for g = 10, ^ approaches inhnity in the limit as g ^ 4, 


Even 


though the q = 7 Potts model exhibits strong hrst-order behavior, correlation length 
is of the order of 50 lattice sites, g = 5 state model has a weak hrst-order transition 
with a correlation length about few thousand in lattice units. In various studies on 
hnite size scaling behavior of the g-state Potts model, it has been shown that g > 10 
state Potts model enters into asymptotic regime for the lattice sizes L > 40 B 1,0 
(for a review see |^). This makes the observation of scaling difficult for g < 10 in 
simulation studies by using computationally feasible lattice sizes B [ 
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higher order corrections become important. However these studies are based on the 
hnite size scaling behavior of the local operators such as the energy cumulants. The 
aim in this work is to observe the characteristic behaviours in the system by studying 
global operators such as huctuations in cluster size which may yield information on 
the scaling form of the operators and hence on the order of the transition. 
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The Hamiltonian of the Potts model |^, [T^ is given by 


(3H = KY. 


( 1 ) 


<*J> 


Here K = J/kT ; where k and T are the Boltzmann constant and the temperature 
respectively, and J is the magnetic interaction between spins Uj and cXj, which can 
take values 1, 2,..., g for the g-state Potts model. One of the basic operators, the 
average cluster size can be defined as 



( 2 ) 


where Ci is the number of spins in the cluster devided by the total number of 


spins. Similar to the energy cumulants, the fluctuations in this quantity at the 
transition temperature take a polynomial form [|l| 


{< > — < C >‘^)max — ^0 + ^i/+ ^ 2 /+ .... 


(3) 


max 


where d is the space dimension. 

In the simulation of 7—state Potts model the measurements are done on twelve 
lattices in the range 12 < L < 64. Since the calculated correlation length ^ of this 
model is about 50 lattice sites [Q , this set of lattices are considered to be an indicative 
set for both very small lattice measurements and for the lattices of the order of the 
correlation length. Measurements are done around the transition temperature and 
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at each temperature 5 x 10® iterations are performed following the thermalization 
runs of 5 X 10^ to 10® iterations. For q = 5 state Potts model, the measurements 
are done on lattice 32 < L < 128. On the largest three lattices 10® iterations were 
performed in order to obtain the same statistical errors. For all simulation works, 
cluster update algorithm |]^ is employed. 


The hrst indication towards the realization of our expectations of the scaling behav¬ 
ior of cluster related operators can be seen in Fig. 1 where cluster size fluctuation 
(FCS) data for g = 7 model is displayed. Using only the lattices with 48 < L, the 
cluster size fluctuations {FCS) data exhibit data collapsing which is a manifestation 
of scaling. Furthermore, as one can see from Fig. 1, the data obtained on lattices 
28 < L < 48, deviate very little around the data collapse curve which means a 
possible correction term in scaling form has very little effect. This is an indication 
that the asymptotic regime is reached for these lattices. Such a data collapse for 
L < 64 is impossible to observe for the energy fluctuations, as it has been studied 
extensively in literature 0. This observation is also supported by the hts done to 
the maxima of the cluster fluctuation data {FCS)max- This quantity is expected to 
be independent of size in the scaling region for a strong hrst-order phase transition. 
In Fig. 2, the {FCS)max values are plotted against The {FCS)max of the 

largest three lattices (L=48, 56 and 64), in most conservative view, are size inde¬ 
pendent. Increasing the number of data points by considering the next two lattices 
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gives the same constant within errors and improves the quality of the ht. The ht 
can be repeated by including smaller size lattices {L < 40) in which case the hrst 
correction term in Eq. (3) must be included (FORMl = Aq + Ai/L*^). When all the 
lattice sizes are included, the results are Aq = 0.0843±0.0003 and Ai = 0.630T0.064 
with = 0.43. The third term in Eq. (3) seems to have no effect on the hrst two 
coefficients and the large error on the coefficient A 2 indicates that this term has no 
room for this set of data. 

In the case of 5—state Potts model, since the correlation length is of the order of 
few thousands of spins, observing asymptotic behavior by using lattices of size in 
the range of L = 32 to 128 sites can not be expected. Nevertheless, if one can see 
the expected hnite size behavior with the correction terms, this can be indicative of 
the order of the transition. Fig. 3 shows {FCS)max versus 1/L'^ for q = 5 Potts 
model. A linear ht in 1/L^ (FORMl) to the largest four lattices (64 < L) gives 
Aq = 0.0429 ± 0.0007 and Ai = 48.1 ± 5.7 with = 0.16. For a ht to the smaller 
lattices, A 2 / term becomes necessary (FORM2 = Aq + Ai/+ A 2 / This 
form hts to the lattices between size 128 down to 56. In this procedure Aq and Ai 
remains unchanged within the error bars. In order to include smaller lattices another 
form (FORM3 = Ao+Ai/L'^+A 2 ea;p(—A 3 L)) has also been tested for all size lattices. 
This form hts to all of the lattices (Fig.3.) while leaving the hrst coeehcient (Aq) 
unchanged. 
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Recently there has been a new pnblication in snpport of the idea that the global 


operators carry more information than the local operators . The hnite size scaling 


algorithm based on bulk and surface renormalization of de Oliveira [18| is tested 


on q-state Potts models in dimensions d = 2 and 3. The operators chosen are very 
sensitive to the order of the transition. They have calculated the magnetic critical 
index and obtained with good accuracy the hrst-order value for even weak hrst-order 
transitions such as for the d = 2, q = 5 Potts model. This is in support of our 


expectations. 


In conclusion, the cluster size fluctuations data for 7-state Potts model on rather 
moderate size lattices hts to the expected behavior of the hrst-order phase transitions 
with scaling exponent being the space dimension. Measurements on lattices as small 
as L = 28 indicate that a very small correction is needed for scaling of the cluster 
huctuation peaks. Furthermore, as it is expected from the large correlation length, 
the q = 5 case can not show such good scaling behavior, but nevertheless the 
expected form for hrst-order phase transition can be observed with very stable hts 
to the data. In the case of energy huctuations, the data collapse, even for the q=7 
state Potts model case, is seen to be impossible and the lattices considered failed 
to give even stable hts. In this respect it can be emphasized that the choice of the 
operators is extremely important, especially in the case of weak hrst-order phase 
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transitions. The resnlts are in favor of previous observations that the operators 
which possess global information on the system can better exhibit the character of 
the phase transition. 
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Figure captions 


Figure 1. FCS data versus {K — Kc) x for q = 7 Potts model. 

Figure 2. {FCS)max plotted as a function of for g = 7 Potts model, 
sizes are 12 < L < 64. 

Figure 3. {FCS)max plotted as a function of for g = 5 Potts model, 
sizes are 32 < L < 128. 
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